JEM 



Article 



The distinctive germinal center phase 

of IgE + B lymphocytes limits their contribution 

to the classical memory response 

Jin-Shu He, 1 Michael Meyer-Hermann, 2 Deng Xiangying, 1 

Lim Yok Zuan, 1 Leigh Ann Jones, 1 Lakshmi Ramakrishna, 1 

Victor C. de Vries, 1 Jayashree Dolpady, 3 Hoi Aina, 1 Sabrina Joseph, 1 

Sriram Narayanan, 1 Sharrada Subramaniam, 1 Manoj Puthia, 1 Glenn Wong, 1 

Huizhong Xiong, 3 Michael Poidinger, 1 Joseph F. Urban, 4 Juan J. Lafaille, 3 

and Maria A. Curotto de Lafaille 1 ' 3 



C 

u 



d 

S3 
o 

s 

U 

o 

W 
0 

73 

53 
sh 

3 

0 

> 

H 



CORRESPONDENCE 

Maria A. Curotto de Lafaille: 

maria_lafaille@immunol.a-star 

.edu.sg 

Abbreviations used: CSR, class 
switch recombination; DZ, dark 
zone; FDC, follicular DC; GC, 
germinal center; GSEA, Gene 
set enrichment analysis; LZ, 
light zone; mLN, mesenteric 
LN; PC, plasma cell; Tfh, 
T follicular helper cell. 



Singapore Immunology Network (SlgN), Agency for Science, Technology and Research (A*STAR), Singapore 138648 
department of Systems Immunology, Helmholtz Centre for Infection Research and Department of Life Sciences, Technische 
Universitat Braunschweig, D-38124 Braunschweig, Germany 

3 The Kimmel Center for Biology and Medicine of the Skirba 1 1 Institute and Department of Pathology, New York University 
School of Medicine, New York, NY 10016 

4 Beltsville Human Nutrition Research Center, Agricultural Research Service, United States Department of Agriculture, 
Beltsville, MD 20705 



The mechanisms involved in the maintenance of memory IgE responses are poorly under- 
stood, and the role played by germinal center (GC) lgE + cells in memory responses is par- 
ticularly unclear. lgE + B cell differentiation is characterized by a transient GC phase, a bias 
toward the plasma cell (PC) fate, and dependence on sequential switching for the produc- 
tion of high-affinity IgE. We show here that lgE + GC B cells are unfit to undergo the 
conventional GC differentiation program due to impaired B cell receptor function and 
increased apoptosis. lgE + GC cells fail to populate the GC light zone and are unable to 
contribute to the memory and long-lived PC compartments. Furthermore, we demonstrate 
that direct and sequential switching are linked to distinct B cell differentiation fates: direct 
switching generates lgE + GC cells, whereas sequential switching gives rise to lgE + PCs. We 
propose a comprehensive model for the generation and memory of IgE responses. 



IgE antibodies are critical mediators of allergic 
reactions (Gould and Sutton, 2008). Cross-linking 
of IgE molecules bound to high affinity FceRI 
receptors on mast cells and basophils leads to 
the rapid release of potent proinflammatory 
molecules (Kinet, 1999; Galli andTsai, 2012). 
In spite of its pathological potential, IgE exhibits 
the lowest serum concentration and the short- 
est half-life of all the antibody isotypes (Vieira 
and Rajewsky 1988; Gould and Sutton, 2008). 
The low frequency of IgE -producing cells makes 
their study particularly challenging. Using 
mouse models of high IgE responses (Katona 
et al, 1988; Curotto de Lafaille et al., 2001), we 
discovered that IgE-producing cells develop via 
a unique differentiation pathway that occurs 
during the germinal center (GC) phase of 
T cell— dependent responses and yet favors the 
production of plasma cells (PCs; Erazo et al., 
2007;Yang et al., 2012). In our early studies a 
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GC IgE + population was not clearly detectable, 
but the IgE antibodies produced were observed 
to have undergone affinity maturation, indi- 
cating a GC history for IgE + PC. We proposed 
at the time that high affinity IgE originated 
from the sequential switching of high affinity 
IgGl cells, and hence we speculated that classi- 
cal IgE + memory cells may be absent in mice 
(Erazo et al., 2007; Curotto de Lafaille and 
Lafaille, 2010). 

Sequential switching of IgG cells to IgE was 
first discovered by the identification of switch(S)7 
region footprints in the S|JL-Se DNA region of 
IgE genes (Matsuoka et al., 1990;Yoshida et al., 
1990;Jabara et al., 1993; Mandler et al., 1993; 
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Zhang et al., 1994; Baskin et al., 1997), but the biological sig- 
nificance of this finding was at that time unknown. Sequential 
switching in mice entails two recombination events, S(x— > 
S7I and S1J1S7I— »Se, that may be either continuous or tem- 
porally separate events. The latter scenario allows for the exis- 
tence of an intermediate IgGl cellular phase in which affinity 
maturation can occur in GCs. Indeed, stimulation of IgGl 
cells in the presence of IL-4 either in vivo or in vitro resulted 
in the production of IgE antibodies (Erazo et al., 2007; 
Wesemann et al., 2012). Importantly, mice deficient in class 
switching to IgGl due to a mutation in the I7I exon (Lorenz 
et al., 1995) were unable to produce high affinity IgE anti- 
bodies (Xiong et al., 2012a,b), indicating that sequential 
switching is essential for the formation of high affinity IgE. 

The recent development of fluorescent reporter mice 
for IgE has facilitated the identification of IgE GC cells (Talay 
et al., 2012;Yang et al., 2012). However, the in vivo phenotype 
and role of IgE GC cells in supporting IgE responses and 
its relationship with the sequential switching process remain 
unclear (Lafaille et al., 2012; Xiong et al., 2012a). 

In the current study, we used a new reporter mouse for 
class switch recombination (CSR) to IgE, improved methods 
to functionally study IgE B cells ex vivo and in vivo, and in 
silico modeling to analyze the origin, functional properties, 
and population dynamics of IgE GC cells and PC. We show 
that IgE GC cells are unfit to undergo the conventional GC 
differentiation program and instead undergo apoptosis at a 
high rate.This "failure to thrive" of IgE GC cells greatly limits 
their contribution to the memory pool and high affinity PC 
compartment. Furthermore, we show that the two types of 
rearrangement to IgE are associated with distinct B cell dif- 
ferentiation fates. Direct S(x-Ss rearrangements generate IgE 
GC cells, whereas sequential switching of IgGl cells gives rise 
to IgE PC. 

RESULTS 

Expression of GFP in CeGFP mice reports all CSR to Ce 

To track immunoglobulin gene CSR to IgE in vivo and 
in vitro, we generated the CeGFP reporter mice. These mice 
carry an IRES-GFP cassette insertion in the 3' untranslated 
region of the membrane-encoding Ce gene that preserves 
native polyadenylation signals (Fig. 1 A and Fig. SI A). In vitro 
stimulation of naive splenocytes from CeGFP mice with ei- 
ther LPS or anti-CD40 in the presence of IL-4 led to the ap- 
pearance of GFP + cells concurrently with production of IgE 
antibodies, whereas LPS or anti-CD40 stimulation alone did 
not (Fig. SI, Band C). 

To assess GFP expression in vivo, we infected CeGFP 
BALB/c mice with the parasite Nippostrongylus brasiliensis, a 
known inducer ofTh2 responses and IgE production (Urban 
et al., 1992). Parasite infection resulted in the expected ap- 
pearance of IgE + GFP + cells and also elicited a population of 
IgGl + GFP + cells in LNs (Fig. 1 B). All GFP + cells expressed 
CD95 (Fig. 1 B). Only IgE + cells expressed the mature IgE 
transcript, whereas the post-switched Ce transcript was ex- 
pressed by both IgE + GFP + cells and IgGl + GFP + cells, but not 
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by IgGl + GFP- cells (Fig. 1 C). It is well estabhshed that CSR 
can occur in both IgH alleles (Rabbitts et al., 1980; Siebenkotten 
et al., 1992; Erazo et al., 2007); thus, the expression of both 
GFP and switched Ce transcript in IgGl + GFP + indicates that 
these cells underwent recombination to Se in the nonproduc- 
tive IgH allele (Fig. 1 D). These data demonstrate that GFP ex- 
pression in the novel IgE reporter CeGFP mice uniquely 
identifies all events of CSR to Ce. 

Marked differences in population dynamics 
of lgE + and IgGl + GC cells 

We next characterized IgGl + and IgE + cells in mesenteric 
LN (111LN) from CeGFP at 10 d after infection with N. brasil- 
iensis. IgE + and IgGl + cells could each be separated into a 
B220 lo CD138/Syndecan-l + PC population that expressed 
Prdml, Irf4, and Xbpl, and a B220+CD138" GC population 
that instead expressed Bel 6 and Aicda (Fig. 1, E and F). Gene 
array analysis revealed very similar gene expression patterns 
within the IgE + and IgGl + (IgGl + GFP + and IgGl+GFP-) 
GC cells (Fig. 1 G). IgE + GC cells expressed lower levels of 
Ig than did IgE + PC cells (Fig. 1 H), whereas IgGl + GC cells 
(both GFP + and GFP - ) expressed higher Ig levels than did 
their IgGl + PC counterparts (Fig. 1 H), as previously reported 
(Erazo et al, 2007; Xiong et al, 2012b;Yang et al, 2012). The 
predominance of the PC in the IgE + population (Erazo et al., 
2007) was evident at 10 d after infection, at which point 
IgE + GFP + PC comprised ~50% of the total IgE + GFP + cells 
in mLN, while IgGl + GC cells greatly outnumbered IgGl + 
PC (Fig. 1 E). 

We next sought to track cells carrying rearrangements 
to Ce during infection with N. brasiliensis, which included 
IgE cells, as well as IgGl cells expressing GFP. Because the mice 
were heterozygous for the CeGFP insertion, GFP + cells com- 
prised approximately half of the total cells carrying Se rear- 
rangements. Due to the low level of IgE surface expression, 
IgE + GFP~ GC cells could not be reliably tracked. IgGl + GFP~ 
and IgGl + GFP + GC cells displayed similar kinetics, albeit 
very different frequency. IgE + GFP + GC expanded quickly up 
to 13 d, similarly to IgGl + GC cells, but declined very rapidly 
thereafter (Fig. 2 A). The similar frequencies of IgE + GFP + 
and IgGl + GFP + GC cells up to 13 d after infection suggests 
that switching to IgE occurred at comparable rates in both 
the productive and nonproductive IgH alleles (Fig. 2 A). As the 
frequency of IgE + GC cells began to decline from day 13 on- 
ward, the IgGl + GFP + cells began to outnumber the IgE + GFP + 
population. During secondary infection, the increase in IgE + 
GC cell frequency was also transient whereas IgGl + GC cell 
numbers were more stable (Fig. 2 A). 

The frequency of IgE + PC (both GFP + and GFP") at 10 d 
after infection was similar to that of IgE + GFP + GC cells, 
but the IgE + PC population decreased rapidly thereafter and 
their decline preceded the loss of IgE + GFP + GC cells (Fig. 2, 
A and B). These data are inconsistent with the previously 
proposed model in which the primary cause of declining 
IgE + GC cell numbers was attributed to be their differentia- 
tion into PC (Yang et al., 2012). After secondary infection 
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Figure 1. CeGFP mice efficiently report all CSR to Ce. (A) Schematic representation of CeGFP Kl integration into the IgHCs locus. (B-H) CeGFP 
BALB/c mice were infected with N. brasiliensis and mLNs were collected 10 d later for analysis by flow cytometry and sorting. (B) Left: frequency of GFP+ 
cells among total lymphocytes; middle: CD95 expression on GFP+ cells and GFP~B220 + cells; right: frequency of lgE + and lgGl + cells within the gated 
GFP+ population. (C) QPCR analysis of Ce switched (left) and mature (right) transcripts in sorted lgE + GFP+, lgG1+GFP+, and lgGl + GFP~B220 + cells. Error 
bars indicate SD of duplicated samples. Data are representative of two independent experiments with n = 3 mice pooled per experiment. (D) Schematic 
representation of mature Ce and switched Ce transcripts driving GFP expression from the productive and nonproductive IgH alleles. (E) Numbers indicate 
frequency of B220'°CD138+ PC and B220+CD138- GC cells among gated lgE+GFP\ lgGl + GFP+, and IgG+GFP- cells. (Fand G) IgE+and lgGl+ LN cells from 
day 10 N. brasiliensis CeGFP BALB/c-infected mice were sorted into B220'°CD138 + PC and B220 + CD138~ GCcell populations. Naive B220+FAS lgG1 IgE 
cells were sorted as a control reference population. Gene expression was analyzed using the Affymetrix Mouse Exon 1.0 ST Array. (F) Graphs showed hy- 
bridization intensities for GC genes Bc/6 and Aicda and for PC genes Prdml, Irf4, and Xbpl within the sorted populations. Error bars indicate SEM of n = 3-4 
samples per group. (G) Heat-map representation of the 8,120 probes differentially expressed in at least one comparison between group IV and groups I, II, 
or III in the GC array. Probes have been z-score normalized. Similar patterns of gene expression were observed across groups I, II, and III, but not in group IV. 
(H) Membrane immunoglobulin levels on gated lgE + GFP + , IgGl + GFP + , and IgG+GFP" GC and PC cells. Data shown are representative of at least five experi- 
ments (B, E, and H), n = 3-5 mice per experiment. 



IgE + and IgGl + PC frequency in LN transiently increased 
(Fig. 2 B). IgE + PC and IgGl + PC cells were detected in the 
BM at 21 d after primary infection and also after secondary 
infection (Fig. 2 C). 



These kinetic analyses identified two distinct phases of the 
primary GC B cell response in N. brasiliensis— infected mice: an 
initial infective phase that is characterized by the rapid gener- 
ation and expansion of IgE + and IgGl + GC cells, followed by 
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Figure 2. Population dynamics of lgE + 
and IgG 1 + cells during infection with 
N. brasiliensis. CsGFP BALB/c mice were 
infected with N. brasiliensis and a group of 
mice were reinfected 40 d later. (A and B) mLN 
were collected on days 8, 10, 13, 14, 21, and 30 
after primary infection, and on days 7 and 14 
after reinfection. (A and B) Kinetics of lgE + GFP + , 
lgG1 + GFP + , and IgGl+GFP- GC cells (A) and 
lgE + GFP + , IgE+GFP- lgGl + GFP + , and IgGUGFP- 
PCs (B) during primary and secondary N. brasil- 
iensis infection. The frequencies of each cell pop- 
ulation are given as proportion of total DAPI~ 
viable cells. (C) Frequency of BM lgE + PC and 
IgG 1 + PC on day 21 of primary infection (pri- 
inf) and on day 14 of reinfection (re-inf) with 
N. brasiliensis. UT, untreated mice, n = 3-6 
mice per time point from two independent 
experiments were analyzed (A-C). Error bars 
indicate SEM. 



a second phase commencing after worm expulsion on day 1 1 
after infection (Camberis et al., 2003), at which time IgE + 
GC steeply decline and IgGl + GC cells persist. 

Functional deficiency in B cell receptor 
expression and signaling in IgE GC cells 

We next hypothesized that a deficiency in BCR expression 
and/ or signaling could underlie the surprisingly rapid decline 
of IgE + GC cells in parasite-infected CeGFP mice. To directly 
compare BCR expression between IgE + and IgGl + cells, we next 
studied these populations in TBmc mice (Curotto de Lafaille 
et al., 2001), in which all B cells are specific for the HA pep- 
tide. This enabled us to use the same peptide antigen to stain 
the BCR in both IgGl + and IgE + cells. BCR levels were quan- 
tified in mLN cells isolated on day 12 after immunization by 
staining with the peptide HA and with antibodies against IgE 
and IgGl. Mutations inV(D)J genes are rare at this time point 
in TBmc mice (Erazo et al., 2007); hence, HA binding was nei- 
ther compromised, nor did it interfere with the binding 
of isotype-specific antibodies (Fig. 3, A and B). Analysis of 
HA binding level (MFI) demonstrated that surface immuno- 
globulin levels in IgE + GC cells were three- to fourfold lower 
than in IgGl + GC cells (Fig. 3, C and D). In contrast, IgE + 
PCs expressed three- to fourfold higher levels of membrane 
immunoglobulin than did IgGl + cells (Fig. 3 D). 

We next investigated the correlation of surface BCR ex- 
pression with transcript levels for total Ig, membrane-bound 
Ig, and secreted Ig (Fig. 3 E).We detected comparable levels 
of mature steady-state transcripts in both IgE + and IgGl + GC 
cells, as well as in both IgE + PC and IgGl + PC, although total 
transcript expression level was > 100-fold higher in PC cells 
compared with GC cells (Fig. 3 F). IgE + GC cells produced 
atypically high levels of secreted IgE transcripts but low levels 
of IgE membrane transcripts (Fig. 3, G and H), as has previ- 
ously also been described for IgE + B220 + cells (Karnowski et al., 
2006). In addition, membrane IgE transcripts were up-regulated 
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in IgE + PC compared with IgE + GC cells, whereas the re- 
verse was the case for IgGl + cells (Fig. 3 H), consistent with 
our earlier flow cytometry data (Fig. 1 H). 

BCR expression and signaling are essential for the sur- 
vival of mature B cells (Lam et al., 1997; Kraus et al., 2004), 
and both BCR/Ag-internalization and BCR signaling are 
thought to be necessary for effective selection and survival 
in the GC (Tarlmton and Smith, 2000; Hou et al., 2006; 
Vascotto et al., 2007; Victora and Nussenzweig, 2012). We 
therefore investigated whether the rapid decline of IgE + GC 
cells in CeGFP mice could be due to impaired BCR signal- 
ing in this population. A previous report identified that IgM + 
GC cells exhibit impaired BCR signaling (except during 
G2/M) due to a high level of phosphatase activity (Khalil 
et al., 2012), and similarly, we observed that switched IgGl + 
and IgE + GC cells from TBmc mice responded poorly to 
BCR cross-linking (Fig. 4 A). Treatment of GC cells with 
increasing concentrations of the phosphatase inhibitor H 2 0 2 
led to a dose-dependent phosphorylation of proximal signaling 
molecules Syk and Blnk (Fig. 4, B and C). Phosphorylation 
of Erkl/2 occurs downstream of BCR and CD40 signahng, 
and high levels of p-Erkl/2 expression were induced even at 
the lowest H 2 0 2 concentration (Fig. 4, B and C). Interest- 
ingly, when subjected to suboptimal stimulation, IgE + GC 
cells exhibited far lower levels of p-Syk and p-Blnk than did 
IgGl + GC cells (Fig. 4, B and C), consistent with a reduced 
number of BCR signaling complexes in IgE + cells. It is un- 
likely that this reflects inherent differences in phosphatase ac- 
tivity because both the SHP-1 and SHIP-1 phosphatases that 
regulate BCR signaling (Cyster and Goodnow, 1995; Pani 
et al, 1995; Liu et al., 1998; Okada et al, 1998) were ex- 
pressed at similar transcriptional levels in IgE + and IgGl + GC 
cells (not depicted). 

In addition to the BCR levels, we quantified the expres- 
sion of other surface receptors important for antigen presen- 
tation and costimulation, such as MHC II, CD19, CD40, 
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Figure 3. Dysregulated expression of mem- 
brane and secreted immunoglobulin in lgE + 
cells. (A-D) TBmc mice were immunized with 
0VA-HA and the mLNs were harvested 12 d later. 
Cells were incubated with 0-1,000 ng/ml of labeled 
HA peptide, along with antibodies. (A) Dot plots 
show the frequencies of lgE + and lgGl + cells 
among the gated lgD-|gM-B220 + CD138- switched 
GC cells and among CD138+B220' 0 PC populations. 
(B) Binding competition between HA peptide and 
antibodies against membrane immunoglobulin was 
examined by flow cytometry. Histograms show 
fluorescence intensity of IgE and IgGl staining in 
GC cells costained with increasing concentrations 
of HA. (C) Histograms show levels of surface- 
bound HA on lgE + and lgGl+ GC cells. (D) MFI of 
bound HA on GC populations (top) and on PC pop- 
ulations (bottom) was determined at different 
peptide concentrations by subtracting the MFI 
value of the HA unstained control samples. 

(E) QPCR primer locations for total IgE, membrane- 
bound IgE, and secreted IgE transcripts (tr.). 

(F-H) CeGFP TBmc mice were immunized with 
0VA-PEP1 and the mLN were harvested 12 d later. 

(F) QPCR analysis of total VDJ transcript levels in 
purified lgE + and lgGl + GC cells and lgE + and 

IgG 1 + PC. (G and H) QPCR analysis of secreted and 
membrane Ce transcripts, and secreted and mem- 
brane C7I transcripts in lgE + and IgG 1 + GC cells 
and in lgE + and lgG1 + PC. Error bars indicate SEM 
of n = 3 mice (D and F-H). *, P < 0.05; **, P < 0.01 ; 
***, P < 0.001. Data are representative of two inde- 
pendent experiments (A-D and F-H). 




GC PC 



GC PC 



GC PC 



GC PC 



CD80, LIGHTR, CD21/CD35, ICOSL, and OX40L. We 
found that expression of the surface C3d receptor CD21/ 
CD35, and of ICOSL and OX40L, were significantly re- 
duced in IgE + GC cells compared with IgGl + GC cells (Fig. 5, 
A and B). 



In sum, we demonstrated marked differences in the regu- 
lation of membrane and secreted Ig expression between IgGl + 
and IgE + cells, and in particular, that IgE + GC cells exhibit 
decreased BCR expression and reduced proximal signaling 
upon phosphatase inhibition. 
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Figure 4. BCR signaling is reduced in lgE + GC cells compared with lgG1 + GC cells upon ex vivo stimulation. (A) Single cell suspensions of ml_N 
from OVA-HA immunized CeGFP TBmc mice were treated with D-PBS or stimulated ex vivo with 10 ug/ml OVA-HA for 5 min. Levels of p-Syk and p-BInk 
expression were determined by flow cytometry within gated populations of lgE + GFP + cells, lgG1 + GFP+, and IgG 1 Ig E GC cells, and B220+CD95~ non-GC 
cells (black line histograms). Gray filled histograms correspond to D-PBS-treated controls. Data are representative of two independent experiments (n = 9 
mice pooled per experiment). (B) Cells from immunized CeGFP TBmc mice were untreated or stimulated with increasing concentrations of H 2 0 2 for 
10 min. Histograms show p-Syk, p-BInk, and p-Erk1/2 expression levels in gated lgE + GFP + GC cells (black line) and in lgG1 + GFP~ GC cells (solid gray). Data 
are representative of nine independent experiments [n = 6-9 mice pooled per experiment). (C) Quantification of p-Syk, p-BInk and p-Erk1/2 levels (MFI) in 
H 2 0 2 -stimulated Ig E + GFP + and lgGl+GFP~ GC cells from four independent experiments (n = 6-9 mice pooled per experiment). Error bars represent SEM. 
*, P < 0.05;**, P<0.01. 



Increased dark zone (DZ) apoptosis and reduced 
light zone (LZ) frequency suggest decreased 
DZ to LZ output of lgE + GC cells 

GCs are anatomically and functionally separated into two 
distinct zones: a DZ, wherein cell proliferation and somatic 
hypermutation occur, and an LZ, where GC cells undergo se- 
lection via interactions with follicular DCs (FDCs) andT follicu- 
lar helper cells (Tfh cells; Kelsoe, 1996;Victora and Nussenzweig, 
2012). The LZ and DZ environments are maintained by dif- 
ferential expression of the chemokines CXCL12 (DZ) and 
CXCL13 (LZ), and GC B cells in these different zones can be 
distinguished by the differential expression of CXCR4 and 
CD86/CD83 (Allen et al., 2004;Victora et al., 2010). 

Intriguingly, analysis of the LZ/DZ distribution of IgE + and 
IgGl + cells demonstrated a threefold reduction in the frequency 
of IgE + cells exhibiting the characteristic CXCR4~CD86 + 
LZ phenotype compared with IgGl + cells (Fig. 6, A and B). 
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Differentially expressed LZ and DZ gene sets were then gen- 
erated by comparison of the LZ and DZ gene expression 
profiles of IgGl + cells (Table S2). GeneGo pathway analysis of 
LZ and DZ genes yielded comparable results to previously 
reported data (Victora et al, 2010, 2012), although we did 
observe a novel increase in mitochondria redox activity in 
DZ cells (Table S2).This pathway was found to be specifically 
up-regulated in GC cells among all B lymphocyte popula- 
tions when compared in Immgen (Heng and Painter, 2008). 
Gene set enrichment analysis (GSEA) of the LZ and DZ 
gene signatures of total IgE + and IgGl + GC cells (Fig. 6 C 
and Table S2) revealed enrichment of the DZ gene signature 
among IgE + cells (Fig. 6 C), which was consistent with our 
earlier flow cytometry data. BrdU incorporation identified a 
higher frequency of proliferating cells in the total IgE + GC 
population than in total IgGl + GC cells, reflecting the low 
LZ component of IgE + GC cells (Fig. 6 D). 
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Figure 5. Decreased expression of costimulatory molecules on lgE + GC cells. CeGFP BALB/c mice were infected with N. brasiliensis and mLN cells 
were harvested and analyzed 10 d later. (A) Flow cytometry analysis of MHC II, CD19, CD40, CD80, LIGHTR, CD21/35, 0X40L, and ICOSL surface expression, 
and assessment of total (surface and intracellular) expression of CD21/35, 0X40L, and ICOSL in gated lgE+GFP+ and IgGl+GFP" GC cells (FMO = Fluores- 
cence Minus One control). (B) Quantification of surface expression levels (MFI) of CD21/35, OX40L, and ICOSL on gated IgE+GFP+and lgG1+GFP- GC cells. 
Error bars indicate SEM of n = 5 mice. ***, P < 0.001. Data are representative of three independent experiments (A and B). 



In addition to their low LZ frequency, IgE + cells did not 
display up-regulation of BCR expression in the LZ relative to 
their counterparts in the DZ, whereas BCR levels among 
IgGl + cells were ^40% higher in LZ cells than in DZ cells 
(Fig. 6 E). These differences in distribution and DZ/LZ BCR 
expression were also evident in GC sections stained with anti- 
bodies against IgE/IgGl and with anti-CD35 to identify 
FDC (Fig. 6 F). As previously reported (Erazo et al., 2007), 
brighter staining of IgGl + cells was observed in the LZ of 
the GC, consistent with higher BCR expression in IgGl + LZ 
cells than in the IgGl + DZ cells. Most IgE + GC cells ap- 
peared to be located in the DZ, reflecting the increased fre- 
quency of these cells in the DZ. 

Because low expression of BCR and related costimula- 
tory molecules may impair the selection and survival of GC 
cells, we next investigated whether IgE + cells undergo apop- 
tosis at a higher rate than do IgGl + cells. Analysis of total 
Caspase activity among these populations (Fig. 6, G and H), 
together with assessment of activated Caspase-3 levels (not 
depicted), consistently indicated a two- to threefold higher 
incidence of apoptosis among IgE + GC cells compared with 
IgGl + GC cells (Fig. 6 H). Unexpectedly, the majority of 
IgE + apoptotic cells (68%) accumulated in the DZ, whereas 
the majority of IgGl + apoptotic cells (65%) were located in 
the LZ (Fig. 6,1 and J). 



The distribution of LZ and DZ B cells is determined by 
GC dynamics of cell proliferation, selection, apoptosis, dif- 
ferentiation, and intra-zonal migration (Victora et al., 2010; 
Meyer-Hermann et al., 2012). To determine the best-fit hy- 
pothesis for the IgE DZ/LZ distribution identified here, we 
next developed a mathematical model of CSR and IgE + GC 
B cell properties using previously defined parameters for GC 
dynamics (including cell migration, division, selection, zonal 
exit, and antibody production; Meyer-Hermann et al., 2006, 
2012; Figge et al, 2008; Garm et al., 2010; Zhang et al, 2013). 
This model also incorporated new parameters for IgM + GC 
population kinetics, probability of CSR to IgE and IgGl, and 
immunoglobulin isotype-specific BCR expression. 

Based on the reduced BCR expression of IgE + GC B cells, 
two non— mutually exclusive scenarios were considered in silico. 
In the first scenario, IgE + GC cells were less competitive than 
IgGl + GC cells for antigen capture and were thus impaired 
in their ability to interact withTfh cells and undergo selection 
in the LZ. Inefficient selection of IgE + GC cells in the LZ was 
tested by reducing the probability of antigen uptake by IgE GC 
cells to <0.3 in the in silico model (based on BCR expression 
level among IgE + cells being only 30% of that observed in 
IgGl + GC, Fig. 3 D). However, this inefficient selection pro- 
cess was not sufficient to reproduce the DZ/LZ ratio of IgE + 
cells observed m the wet laboratory (Fig. 6 K). In contrast, in the 
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Figure 6. Altered LZ/DZ distribution and 
increased apoptosis of lgE + GC cells com- 
pared with IgGl + GC cells. Single cell suspen- 
sions were prepared from the ml_l\l of CeGFP 
BALB/c mice on day 10 of infection with 
N. brasiliensis. (A) Dot plots show the DZ and 
LZ phenotypes among lgE + GFP + , lgGl + GFP + , 
and IgGI+GFP GC cells. CXCR4 and CD86 were 
used as markers to discriminate GC LZ 
(CD86 hi CXCR4'°) and DZ (CD86 lo CXCR4 hi ) cells. 
(B) Quantification of DZ to LZ cell ratio among 
the gated populations. Error bars indicate SEM 
of n = 6 mice. (C) GSEA analysis showing the 
enrichment of gene signatures for DZ genes in 
lgE + cells and for LZ genes in IgGl + cells. Gene 
sets were derived from the DEG of the LZ/DZ 
array. All nominal P-values and FDR rates were 
<0.001. (D) Kinetics of BrdU incorporation in 
total (top) and DZ (bottom) lgE+GFP+ and 
IgGI+GFP - GC cells as assessed by flow cytom- 
etry. Error bars indicate SEM of n = 3 mice. 
(E) Ratio of LZ/DZ BCR expression levels (MR) on 
gated LZ (CD86 h CXCR4'°) and DZ (CD86'°- 
CXCR4 hi ) lgE+GFP+ GC cells and lgG1+GC (GFP+ 
and GFP~) cells. Error bars indicate SEM 
of n = 6 mice. (F) Immunohistology of frozen 
mLN sections from day 1 1 after infection with 
N. brasiliensisto show the distribution of lgE + 
and lgG1 + cells (green) in the LZ (CD35 + , red) 
and DZ of the GC. Bar, 50 urn. (G-J) Apoptosis 
assay performed using the CaspGLOW active 
Caspase Staining kit. (G) Flow cytometry analy- 
sis showing the frequency of apoptotic cells 
(CaspGL0W + ) among the gated populations. 

(H) Quantification of cell apoptosis (Casp- 
GL0W+). Error bars indicate SEM of n = 5 mice. 

(I) Apoptotic cell distribution in the CXCR4 hi and 
CXCR4' 0 compartments of the gated popula- 
tions. (J) Percentage of apoptotic DZ cells among 
total apoptotic cells (DZ plus LZ) in lgE + GFP + and 
lgG1 + GFP- GC cells. Error bars indicate SEM of 

n = 5 mice. Data are representative of at least five 
independent experiments (A, B, and E), two inde- 
pendent experiments (D and F), and three inde- 
pendent experiments (G-J). * P < 0.05; **, 
P < 0.01 ; ***, P < 0.001. (K and L) DZ to LZ ratio 
in silico. (K) DZ to LZ ratio when selection of IgE 
GC cells in the LZ was reduced to 30°/o. (L) DZ to 
LZ ratio when IgE DZ cells failed to up-regulate 
CXCR5 and migrate to the LZ. (M) Frequency of 
total apoptotic lgE+ and lgG1+ GC cells (left), 
apoptosis in the DZ (middle), and apoptosis in 
the LZ (right) resulting from impaired DZ to LZ 
output in silico. Data represent mean and SD of 
30 simulations (K-M). 



second scenario where IgE + DZ cells failed to respond to LZ- 
derived CXCL13 and were thus retained in the DZ, the experi- 
mentally observed DZ/LZ distribution of IgE + cells was 
successfully replicated by the in silico model (Fig. 6 L). In accor- 
dance with the experimental data, the reduced DZ to LZ output 
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of IgE + cells in silico was associated with increased rates of apop- 
tosis among this population and similar percentage of apoptotic 
IgE + cells in LZ and DZ (Fig. 6 M).As the DZ/LZ ratio of IgE + 
cells is ^6 (Fig. 6 B), this percentage distribution is consistent 
with most apoptotic IgE + cells residing in the DZ. 
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Figure 7. Direct switching origin of lgE + GC 
cells: molecular and in silico analysis. (A) Sche- 
matic representation of S7I footprint analysis. The 
switched Sp.-Se DNA fragments were PCR amplified 
from purified cells that contained rearrangements 
to Se (lgE+ and lgGl+GFP+) and the PCR products 
were subsequently cloned. Clones containing Su,-Se 
switch regions were identified by QPCR analysis of 
(ji-enhancer and Se sequences. The presence of S7I 
sequences was determined by QPCR or DNA se- 
quencing. (B) S-/1 DNA repeat "footprints" in the 
Sjjl-Se fragment of lgE + cells. S7U clones among 
S|x-Ss DNA clones in lgE + cells sorted from mice 
after primary and secondary infection/immuniza- 
tion. Mice were subjected to secondary infection at 
30-32 d after primary infection or secondary im- 
munization at day 35 after primary immunization. 
lgE + GC cells were sorted from pooled mLN and 
spleen, and lgE + PCs were sorted from pooled mLN 
and spleen or BM. n = 3-7 mice per group total 
pooled from two to three independent experiments. 
(C) Direct switching model in silico. The graph 
shows the population kinetics of lgE + GC cells when 
they were originated from lgM + but not lgG1 + GC 
precursor cells. (D) Sequential switching model in 
silico. The graph shows the population kinetics of 
lgE + GC cells if they were originated from lgG1 + but 
not lgM + GC precursor cells. Data represent mean 
and SD of 30 simulations (C and D). 



Both the experimental and in silico analyses of IgE + GC 
B cells reported here support a model in which reduced DZ 
to LZ output of IgE + cells drives the distinct DZ/LZ distri- 
bution of this population. Reduced expression of the BCR 
and other costimulatory molecules in IgE + GC cells may si- 
multaneously contribute to the impaired selection of these 
cells in the LZ. 

Direct versus sequential switching origin 
of lgE + GC cells and lgE + PC 

We previously demonstrated that sequential switching is the 
main mechanism to generate high affinity IgE antibodies 
(Erazo et al., 2007; Xiong et al., 2012b), so we next investi- 
gated whether direct switching and sequential switching had 
different roles to play in the generation of IgE + GC cells and 
IgE + PC. Genetic evidence of S|xS~yl — >Se sequential switch- 
ing can be determined by the presence of S~yl DNA rem- 
nants in the S|x-Se region of IgE + cells (Xiong et al., 2012b; 
Fig. 7 A). Analysis of S|x-Se DNA fragments in GC cells 
obtained during primary and secondary responses to N. brasil- 
iensis infection (BALB/c mice) or immunization (TBmc mice) 
demonstrated an almost complete absence of S7I remnants 
(Fig. 7 B). In contrast, S~yl + clones were consistently detected 
among S|JL-Se DNA clones from IgE + PC harvested during 
primary and secondary responses (Fig. 7 B).The highest fre- 
quency of S7I remnants (^50%) was found in BM IgE + PC 
after secondary immunization (Fig. 7 B). It should be noted 
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here that only a fraction of S|xS"yl — >Se junctions retain S7I 
DNA (Xiong et al., 2012b); hence, the percentages indicated 
above underestimates the true extent of sequential switch- 
ing. These results indicated that IgE + GC cells are generated 
by direct switching of IgM + cells, whereas sequential switch- 
ing from IgGl + cells gives rise to IgE + PC, even during pri- 
mary immunization. 

We next investigated the implications of these data by 
testing different models of CSR under distinct sets of GC 
kinetic conditions in silico.The population dynamics required 
to induce a steep decline in IgE + GC cells could only be 
achieved when the IgE + GC cells originated from IgM + pre- 
cursors (Fig. 7 C) but not IgGl + precursors (Fig. 7 D). In this 
model, the IgM + cells switched to either IgGl or IgE with 
assigned probabilities of. Pmg 

= 0.04 and 

Pme ~ 0.01, respec- 
tively. IgGl + GC cells were assumed to sequentially switch 
to IgE with probability p GE . In simulations with p GE > 0, the 
population kinetics of IgE + GC B cells matched the popula- 
tion kinetics of the IgGl + cells, even when IgE + B cells failed 
selection at a rate of 100%. Thus, the mathematical analyses 
predicted that IgE + GC B cells are derived from IgM + B cells, 
and that IgGl + GC B cells do not switch to IgE before exit- 
ing the GC. Collectively, both the molecular and in silico 
analyses support a direct switch origin for IgE + GC cells. 

We then determined whether IgE + and IgGl + GC and 
PC cells displayed different rates of somatic hypermutation 
and affinity maturation. Analysis of the heavy chain VDJ gene 
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Figure 8. lgE + GC cells are capable of acquiring 
affinity enhancing maturations. (A) Analysis of VDJ 
nucleotide mutations in sorted lgE + and lgG1 + cells 
from TBmc mice after primary (day 12) and second- 
ary (day 7) immunization with 0VA-PEP1. Mice were 
subjected to secondary immunization 33 d after pri- 
mary immunization. The pie charts show the propor- 
tion of sequences carrying the indicated number 
of nucleotide mutations/sequence {n = 45-59 se- 
quences per group). The mean number of mutations 
is shown in the center of each pie chart. (B) Analysis 
of VDJ amino acid mutations in sorted lgE + and lgGl + 
cells from TBmc mice after primary (day 12) and sec- 
ondary (day 7) immunization with 0VA-PEP1. The 
bar graph shows the mean number of amino acid 
mutations per sequence. Error bars represent SEM 
of n = 45-59 sequences per group. * P < 0.05; 
***, P < 0.001. (C) The table shows the frequency of 
CDR3 high affinity mutations (in red). 
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mutations in IgE + and IgGl + GC cells indicated higher fre- 
quencies of nucleotide and amino acid mutations during pri- 
mary immunizations in IgGl + GC cells compared with IgE + 
GC cells (Fig. 8, A and B). IgE + and IgGl + GC cells carried 
slightly higher numbers of nucleotide and amino acid muta- 
tions per sequence than did PC after primary immunization, 
whereas the reverse was true for IgGl + cells during secondary 
immunization (Fig. 8, A and B).The frequency of high affinity 
amino acid mutations was increased in IgGl + PC compared 
with IgE + PC (Fig. 8 C), consistent with previous observa- 
tions (Erazo et al., 2007;Yang et al., 2012). Our results dem- 
onstrate that IgE + GCs are not intermediates in the generation 
of high affinity IgE + PC through sequential switching, and 
that IgE + GC cells are not hampered in their capacity to ac- 
quire affinity enhancing mutations. 

lgE"B220 + switched B cells and lgE + PC 
are the main contributors to IgE memory 

The impaired GC phase of the IgE response led us to hy- 
pothesize that conventional IgE memory cells might not exist 
(Erazo et al., 2007; Xiong et al., 2012a,b), a concept challenged 
in recent publications (Lafaille et al., 2012;Talay et al., 2012). 
Because markers for memory IgE + cells could potentially differ 
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from markers of other memory B cell populations (Kurosaki 
et al., 2010), we next established a marker-independent method 
for determining whether IgE + B220 + CD138~ cells (contain- 
ing both GC and putative memory cells) could contribute to 
a secondary antibody response. For this purpose we purified 
HA-specific, IgE-depleted Ig]VrigD-B220 + CD138- cells 
and IgE-containing total IgM"IgD^B220 + CD138~ cells from 
TBmc mice immunized 3 wk earlier (Fig. 9 A; and Fig. S2, 
A and B). The purified B cells were transferred to BALB/c mice 
together with naive OVA-specific T cells, and the mice were 
subsequently immunized with OVA-HA. Two weeks later, 
serum levels of HA-specific IgE were assayed and found to be 
comparable between the mouse groups that received switched 
B220 + B cells irrespective of whether the transplanted cells 
included IgE + B cells (Fig. 9 B).The two groups of mice also 
exhibited similar levels of HA-specific IgE transcripts in LN, 
spleen, and BM (Fig. 9 C). Under these conditions, the con- 
tribution of endogenous lymphocytes to the HA-specific IgE 
and IgGl responses was undetectable (Fig. 9, B and C). These 
results indicate that the contribution of IgE + B220 + CD138~ 
cells to the IgE recall response is negligible. 

We next investigated whether IgE-depleted PC from 
spleen and LN were capable of generating an IgE response 
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Figure 9. IgE recall responses depend 
primarily on de novo sequential switching 
rather than on memory lgE + B cells. 

(A) Schematic of transfer experiments in B 
and C. TBmc mice were immunized with 
OVA-HA in alum and switched B cells were 
isolated from pooled spleen and mLN on day 21 
after immunization. 2 x 10 6 switched B cells 
(B220 + CD138-|gM-|gDi or IgE-depleted 
switched B cells were injected i.v. into irradi- 
ated BALB/c mice together with 1 x 10 5 naive 
OVA-specific CD4+ T cells. Recipient mice and 
the control mice that received CD4 + T cells 
only (or no donor cells at all) were immu- 
nized with OVA-HA in alum. (B) Antigen- 
specific serum IgE and IgGl antibody levels in 
recipient mice as assessed by ELISA 2 wk after 
transfer. Error bars indicate SEM. (C) QPCR 
measurement of mature Ce and C7I tran- 
scripts in mLN, spleen, and BM of recipient 
mice 2 wk after cell transfer. Each symbol 
represents an individual mouse; horizontal 
lines indicate the mean, n = 4-7 mice per group 
from three independent experiments (B and C). 



comparable to their CD138~ counterparts. Purified IgE~ PC 
or IgE + PC cells from immunized TBmc mice were trans- 
ferred into BALB/c mice that received no further treatment 
(Fig. 10 A; and Fig. S2, C and D). Two weeks later, the mice that 
had been injected with purified IgE + PC displayed significantly 
higher levels of serum HA-specific IgE antibodies and HA- 
specific IgE transcripts in spleen and BM than mice that had 
received IgE - PC (Fig. 10, B and C). IgE antibodies were still 
detected in serum up to 60 d after transfer (Fig. 10 D). IgE - PC 
transferred production of HA-specific IgGl antibodies but 
not production of HA-specific IgE antibodies (Fig. 10 B), con- 
firming that IgGl + PCs are terminally differentiated and un- 
able to undergo CSR. Although donor IgE + PC contributed 
markedly less to HA-specific serum Ig than did donor IgGl + 
cells, these data clearly demonstrate that IgE + PCs mediate the 
long-term production of serum IgE. Our results are consis- 
tent with the concept that classical memory IgE cells are 
mostly absent in mice, and that IgE memory rehes predomi- 
nantly on de novo switching from non-IgE + cells and on the 
maintenance of long-lived IgE + PC. 

DISCUSSION 

In this study, we investigated CSR to IgE in vivo using a new 
reporter mouse strain that marks all cells carrying rearrange- 
ments to Se. We described novel findings 011 the develop- 
mental origin, phenotype, and function of IgE + GC cells that 
explain their rapid decline and failure to contribute to the 



memory of IgE responses, and the specific role played by 
IgGl + cells in the ultimate generation of IgE + PC. 

We observed striking differences between the population 
kinetics of IgE + and IgGl + GC B cells during a primary im- 
mune response, similar to previous findings from Yang et al. 
(2012). However, although these authors proposed that the 
rapid decline in IgE + GC cell frequency was likely due to their 
differentiation to PC, we now provide evidence for an alterna- 
tive model in which IgE + GC cells are unfit to engage in the 
dynamics of the mature GC and are rapidly depleted by high 
rates of apoptosis.We propose that this phenotype is primarily 
conferred by impaired BCR expression and function in IgE + 
GC cells. In addition, we identified striking differences in the 
switch origin of IgE + GC cells and IgE + PC that further ques- 
tion the ability of IgE + GC cells to give rise to IgE + PC. 

Our analysis of S|x-Ss switch junctions revealed that IgE + 
GC cells are the product of direct class switching, with virtu- 
ally no S7I remnants evident in the S|x-Se junctions after 
either primary or secondary responses. The IgE + GC cells 
present after secondary immunization exhibited more exten- 
sive VDJ mutations than after primary immunization. One 
possibility to explain these findings is that IgM + memory cells 
(Dogan et al, 2009; Pape et al, 2011), rather than IgM naive 
cells, generate IgE + GC cells in secondary immunization. The 
IgM memory would respond faster than naive cells to the 
secondary challenge. In addition (and non— mutually exclu- 
sively) the newly generated IgE + GC cells may rapidly acquire 
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Figure 10. Evidences of BM-homing 
lgE + PCs. (A) Schematic of transfer experi- 
ments in B-D. TBmc mice were immunized 
with 0VA-HA in alum and 10 5 lgE+ PC 
(B220'°CD138 + ) or IgE switched PC 
(lgM-B220'°CD138 + ) were sorted from 
pooled spleen and mLN on day 21 after im- 
munization and injected i.v. into irradiated 
recipient BALB/c mice. (B) HA-specific IgE and 
lgG1 serum antibody levels in recipient mice 
2 wk after cell transfer as determined by 
ELISA. Error bars indicate SEM. (C) QPCR mea- 
surement of mature Ce and C7I transcripts 
in mLN, spleen, and BM from recipient mice 
at 2 wk after transfer. Each symbol represents 
an individual mouse; horizontal lines indicate 
the mean value. (D) Serum HA-specific IgE 
levels at 60 d after transfer as determined by 
ELISA. Serum samples were subjected to two- 
fold predilution. Error bars indicate SEM. Data 
are representative of two independent experi- 
ments (n = 4-5 mice per group per experi- 
ment; B-D). **, P < 0.01 ; ***, P < 0.001. 



So 

< 9. 



IgE 4 
igE- 



iM-m- 



PC 
PC 



mutations due to enhanced T cell help and the presence of 
preexisting GC structures. Alternatively, it is also possible that 
a few IgE + GC cells survive after the primary immunization 
and can transiently expand during a recall response. 

In contrast to IgE + GC cells, S7I remnants were identi- 
fied in 20—30% of the S|x-Se junctions in IgE + PC in spleen 
and LN, and in up to 50% of IgE + PC in BM, indicating a se- 
quential switching origin for these populations. Because S7I 
sequences may be lost during SixS^l-Se switching, the fre- 
quency of sequential switching is always partially underesti- 
mated by the quantification of S"yl remnants. Sequentially 
switched PC cells can never be generated by directly switched 
GC cells because the entire S-yl DNA is deleted upon direct 
SjJL-Se switching. A substantial fraction of the IgE-producing 
PC cells must therefore have distinct origins from those of 
IgE + GC cells. Based on the population kinetics of these sub- 
sets, similar conclusions on the precursor origins of IgE + GC 
cells and IgE + PC were also reached via our mathematical 
modeling approaches. 

We report here that IgE + GC cells express three- to four- 
fold lower levels of BCR than do IgGl + GC cells, as well as 
exhibiting a comparatively reduced BCR signaling response 
to H 2 0 2 .We propose that these deficiencies contribute to the 
high levels of apoptosis observed among IgE + GC cells. BCR- 
mediated antigen uptake is essential for antigen presentation 
to Tfh cells and subsequent clonal selection in the GC (Victora 
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and Nussenzweig, 2012). BCR signaling is also believed to be 
necessary during the GC reaction, although it is not yet clear 
how and when this occurs. Both processes (BCR-mediated 
antigen uptake and BCR signaling) share a requirement for 
signaling molecules such as Syk (Hou et al., 2006;Vascotto et al., 
2007). It has been suggested that BCR signaling-dependent 
checkpoints take place in the G2/M phase during each repli- 
cation cycle (Khalil et al., 2012). Because GC cells in G2/M 
are in the DZ (Victora et al., 2010), BCR signaling would 
occur in the DZ.Thus, impaired signaling in IgE + DZ cells 
may drive the higher incidence of apoptosis that we observed 
among these cells. 

DZ cells may also require BCR function (antigen capture, 
signaling, or antigen presentation) to successfully up-regulate 
CXCR5 and transition into the LZ.We observed that strikingly 
few IgE + cells in the GC exhibited a characteristic LZ pheno- 
type. The DZ and LZ composition of GCs is determined by 
rates of B cell proliferation, selection, differentiation, and inter- 
zonal migration (Meyer-Hermann et al., 2006, 2012; Figge 
et al., 2008; Garni et al., 2010;Victora et al., 2010; Zhang 
et al., 2013). LZ cells are selected to differentiate into memory 
or PCs, or to recycle to the DZ, whereas nonselected cells 
undergo apoptosis. In the DZ these cells instead mutate and 
undergo rapid clonal expansion. A higher proportion of DZ 
than LZ cells engages in intra-zonal migration, producing a 
net migration of cells from DZ to the LZ (Victora et al., 2010). 
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Using mathematical analyses, we tested whether the low fre- 
quency of IgE + LZ cells could be accounted for by impaired 
LZ selection into the recycling pool, or by decreased DZ to 
LZ migration. We found that the best-fit model for the ob- 
served DZ/LZ distribution of IgE + GC cells was generated 
by a decreased DZ to LZ output. It has previously been re- 
ported that low-affinity GC B cells are retained in the LZ 
during competition with high affinity B cells (Victora et al., 
2010). The phenotype of IgE + GC cells therefore appears to 
result from an intrinsic defect in fitness/survival rather than a 
lack of affinity maturation, or differentiation into PC. 

The high apoptosis of IgE + GC cells may partially explain 
the dominance of the PC component in IgE responses. How- 
ever, as previously suggested, an increased tendency of IgE + 
cells to differentiate into PC may also play a role, and we have 
shown that sequential switching from IgGl + cells is coupled 
to IgE + PC differentiation. IgE + PC comprise ^20-40% of the 
total PC population in spleen and LN, similar to IgGl + PC. 
Because class switching to IgE occurs with much lower proba- 
bility than class switching to IgGl (Siebenkotten et al., 1992; 
Hackney et al., 2009), the high frequency of IgE + PC suggests 
a higher rate of differentiation to PC than IgGl + cells. 

The observation in a recent report that the steep decline 
in IgE + GC cells was not prevented that by Bcl2 transgene 
overexpression was considered as proof that IgE + GC cell num- 
bers decrease due to differentiation into PC (Yang et al.,2012). 
However, earlier investigations in Bcl2 transgenic mice re- 
ported that despite exhibiting normal GC size, these animals 
also display restricted proliferation of GC cells in addition to 
impaired GC apoptosis (Smith et al., 1994, 2000). Further- 
more, Bcl2 overexpression led to an increase in the number of 
memory B cells and PC cells with low antigen affinity (Smith 
et al., 1994, 2000). It was concluded that the Bcl2 overexpres- 
sion rescued low affinity GC cells into the memory and PC 
compartments. These complex effects of Bcl2 overexpression 
thus make it difficult to interpret its role in IgE + GC cells. In 
the current report, we demonstrate that IgE + GC cells undergo 
affinity maturation, but their numbers are decreased in the 
LZ where Tfh-mediated selection into memory cells and high 
affinity PC occurs. Thus, IgE + GC cells in the LZ may not 
reach the differentiation stage at which they might benefit 
from transgenic Bcl2 expression. 

To determine whether IgE + GC cells contribute to the 
memory pool and long-lived PC compartments, we per- 
formed adoptive transfer experiments using B220 + switched 
B cells either containing or lacking IgE + cells. We determined 
that IgE + B220 + cells did not significantly contribute to second- 
ary IgE antibody responses in transplanted animals. Rather, 
IgE antibodies were generated by the sequential switching of 
IgGl + cells. It was previously reported that the production of 
IgE, but not IgGl, during secondary immune responses re- 
quires IL-4, suggesting that de novo switching to IgE is neces- 
sary for recall responses (Finkelman et al., 1988). Contrasting 
with these results, Talay et al. (2012) used an IgE-IRES-GFP 
reporter mouse to determine that IgE + memory cells, and not 
other B cells, mediated secondary IgE responses. However, 



the model used by these authors included the addition of a 
human membrane exon and the replacement of key polyaden- 
ylation sites in the IgE locus that may have modified antibody 
expression and/or altered the biology of IgE-producing cells 
in these animals (Lafaille et al, 2012). 

The existence of long-lived IgE + PC has been controver- 
sial (Achatz-Straussberger et al, 2008; Luger et al, 2009;Yang 
et al, 2012), and our current findings now demonstrate that 
IgE + PC contribute to long-term production of serum IgE, 
albeit less efficiently than do IgGl + PC. This may be due to 
increased apoptosis or inefficient BM homing of IgE + PC 
(Achatz-Straussberger et al., 2008), the lower half-life of serum 
IgE (Vieira and Rajewsky, 1988), and cytophilic removal of 
serum IgE through binding to FceRI and FceRII/CD23 
(Gould and Sutton, 2008). 

On average, IgE + PC exhibited fewer high affinity amino 
acid replacements than did IgGl + PC, which was in agree- 
ment with a delayed affinity maturation of the IgE response 
compared with the IgGl response (Erazo et al., 2007;Yang et al, 
2012). If most IgE + PC cells are generated from IgGl + cells, 
it may be that the IgGl + cells that become IgGl + PC are differ- 
ent (of higher affinity or different differentiation stage) than 
the IgGl + cells that become IgE + PC. 

In sum, the data presented here on IgE GC cell pheno- 
type and function, and on the distinct roles of direct and se- 
quential switching in IgE + GC and IgE + PC differentiation, 
can be integrated with previous findings to construct a model 
ofT cell— dependent IgE responses. In the initial phase of a 
response, direct S|x-Se switching from IgM + cells generates 
a transient wave of IgE + GC cells. A population of IgE + PC is 
also generated at least in part through sequential switching 
of IgGl + cells. As total antigen levels decrease and antigen be- 
comes restricted to FDC depots, the generation of IgE + GC 
cells is halted. Reduced BCR function renders IgE + GC cells 
unfit to engage in the GC dynamics, resulting in higher rates 
of apoptosis and a rapid decline of the entire IgE + GC popu- 
lation. Consequently, IgE + GC cells are unable to contribute 
to the memory and high affinity PC compartments. In contrast, 
IgGl + cells can successfully engage in continuous cycles of 
proliferation/mutation and selection in the GC. Interactions 
of IgGl + GC cells or IgGl + memory cells with IL-4— producing 
T helper cells lead to sequential switching to IgE and rapid 
formation of IgE + PC. This model offers an explanation as to 
why the high affinity IgE response depends on sequential 
switching from IgGl + B cells despite the earlier generation of 
IgE + GC cells. It is tempting to speculate that these features 
of the IgE response may result from evolutionary pressure to 
restrict the development of IgE memory cells and yet main- 
tain the capacity to produce limited amounts of high-affinity 
IgE antibodies. 

MATERIALS AND METHODS 

Mice. CeGFP mice were generated by gene targeting in mouse embry- 
onic stem (ES) cells. The DNA targeting vector contained an IRES-GFP 
cassette inserted into the 3' UTR region of the Cfi gene, downstream of 
the M2 membrane exon and upstream of the polyadenylation sequences 
(Fig. SI A). The targeting construct also contained a neo'-hxP cassette inserted 
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in the intron between the CH4 and Ml exons, and an HSV-TK (Herpes 
Simplex Virus thymidine kinase) cassette at the 5' end of the construct. The 
linearized vector was transfected into a clone of the 129-derived Jl ES cell 
line. The ES cells were then selected in medium containing G418 and gany- 
clovir for positive and negative selection, respectively. Colonies carrying 
homologous integrations were first identified by PGR using primers 5'-GAC- 
ACCCCAAGAGAAGGACA-3' (INT5F3) and 5'-GCCTGAAGAAC- 
GAGATGAGG-3 ' (INT5R3), and 5'-CCCACCATCAACAGGGCTA-3' 
(UpM-R), and integration was confirmed by Southern blotting. Retention 
of the IRES-GFP cassette downstream of the membrane exons was con- 
firmed by PGR using the primers 5'-CTTGTGTAGCGCCAAGTGC-3' 
(INT5R) and 5'-TTGCATTCCTTTGGCGAG-3' (I RES R2). The neo T -hxP 
cassette was then removed by Cre-mediated recombination before injection 
of the ES cells into blastocysts. Deletion of the neo'-loxP cassette was con- 
firmed by PGR using primers 5 '-GTATATGTGTTCCCACCACCAGA-3 ' 
(INT5F4) and 5'-CTGTGGTTTCCAAATGTGTCAG-3' (INT5R4),and 
by the loss of resistance to G418. Mice carrying the IRES-GFP KI in the 
Ce gene, before now referred as CeGFP, were backcrossed to BALB/ c mice 
andTBmc mice (Curotto de Lafaille et al., 2001). All the mice used in the 
current study were CeGFP heterozygous. 

BALB/c mice were purchased from The Jackson Laboratory. Mice were 
bred and housed in the specific pathogen-free animal facility of the Biologi- 
cal Research Center (BRC) A*STAR, Singapore. All animal procedures were 
approved by the BRC/A*STAR Institutional Animal Care and Use Com- 
mittee. 6— 8-wk-old mice were used in all experiments. 

Isolation of mouse naive B cells and in vitro induction of IgE pro- 
duction. Naive B cells were isolated from the spleens and mLN of untreated 
CeGFP mice, using anti-CD43 MicroBeads (Miltenyi Biotec) and magnetic cell 
sorting. A total of 0.25 x lOVml of purified cells were stimulated with 20 (j.g/ml 
LPS (Sigma Aldrich) and 10 ng/ml IL-4 {R&D Systems) or with 3 ug/ml anti- 
CD40 antibodies (eBioscience) and 5 ng/ml IL-4 for 4—5 d in vitro. 

Immunizations and infections. CeGFP TBmc mice and TBmc mice 
were immunized by i.p. injection of 100 (j.g OVA-HA (chicken ovalbumin 
cross-linked to HA peptide YPYDVPDYASLRS) or 100 ug OVA-PEP1 
(chicken ovalbumin cross-linked to the low affinity PEP1 peptide YPYDVP- 
DFASLRS) in alum (Erazo et al., 2007). CeGFP BALB/c and WT BALB/c 
mice were infected with 500 L3 larvae of N. brasiliensis via subcutaneous 
route as previously described (Camberis et al., 2003). 

Microarray analysis. NA was isolated from sorted B cell subsets using 
Arcturus Picopure (Life Technologies) RNA Purification kits. RNA quality 
was assessed by Bioanalyzer (Agilent Technologies). SPIA cDNA were pre- 
pared from 2 ng of total RNA using the Ovation Pico WTA system (v2; 
NuGEN) kit protocol before fragmentation and labeling according to the 
Encore Biotin kit (NuGEN) protocol. 

Sample hybridization was performed on MoGene 1.0 ST arrays (Affyme- 
trix; accession no. GSE49033). Data were normalized with RMA and differ- 
ential gene analysis was performed using limma. Both of these analyses and 
the X 2 test f° r set overlap significance were performed in R (v2.15.2; Smyth, 
2004). Expression heat-maps were generated with the Multi Experiment 
Viewer (Saeed et al., 2003). GSEA was performed with the GSEA software 
from the BROAD institute using default parameters and gene set-based per- 
mutation tests (Mootha et al., 2003; Subramanian et al., 2005). All data han- 
dling and analysis was facilitated by Pipeline Pilot (Accelrys). 

Flow cytometry analysis and cell sorting. Staining of IgE + cells was 
preceded by a brief acid treatment to remove FceRII/CD23-bound IgE as 
described previously (Erazo et al., 2007). Staining with the anti-IgE antibody 
R1E4 (Keegan et al., 1991) did not require acid treatment because R1E4 
does not recognize IgE bound to CD23 or FceRI (Mandler et al., 1993). 
All staining of IgE + cells for use in functional analyses (signaling and apopto- 
sis assays or adoptive transfer experiments) was performed using the anti-IgE 
antibody R1E4. 
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Single cell suspensions were incubated with 10 ug/tnl anti-CD 16/32 for 
15 min at 4°C before surface labeling with antibody cocktails in staining 
buffer (2% FBS, 4 mM EDTA, and 0.1% NaN 3 ) for 30 min at 4°C. Intracel- 
lular staining was conducted using Cytofix/Cytoperm kits (BD). Cells were 
analyzed using LSR II 4/5-laser flow cytometer (BD) or sorted using a FAC- 
SAria II 4/5-laser cell sorter (BD). Purified cells were used for signaling 
assays ex vivo, for expression analysis using microarrays or QPCR, or for 
adoptive transfer into recipient mice. 

BrdU proliferation assay. Mice were infected with N. brasiliensis for 10 d 
and subjected to i.p. injection with 3 mg BrdU. After 4-8 h, the mLN were 
collected and BrdU incorporation was analyzed using BrdU Flow kits (BD). 

Apoptosis assays. 10 6 freshly isolated mLN cells from infected or immu- 
nized mice were resuspended in 1 ml RPMI-based complete medium, and 
then immediately mixed with Red-VAD-FMK (CaspGlow; BioVision) and 
incubated for 45 min at 37°C in 5% C0 2 . Cells were then stained with anti- 
bodies against B220, CD138, IgE, and IgGl, and then analyzed in a 5-laser 
LSR II apparatus. For the active Caspase3 assay, freshly isolated mLN cells 
from infected or immunized mice were surface-labeled with antibodies 
against B220, CD138, IgE, and IgGl and then fixed, permeabilized, and 
stained for intracellular active Caspase3 and GFR Cells were analyzed by flow 
cytometry using a LSR II apparatus. 

Quantification of B cell receptor levels in IgE and IgGl cell. Single- 
cell suspensions were prepared from the mLN of OVA-HA— immunized 
CeGFP TBmc mice and then incubated with Alexa Fluor 647— anti-IgE 
(R1E4), PerCP-eFluro710-anti-IgGl, PE-anti-CD138, APC-eFluro780- 
anti-B220, eFluro650-anti-IgM, eFluro450-anti-IgD, and biotin-labeled HA 
peptide (0—1,000 ng/ml) at 4°C for 30 min, followed by addition of strepta- 
vidin-PE.Cy7.The samples were then analyzed using a LSR II apparatus. 

Ex vivo stimulation and analysis of phospho- signaling proteins. Sin- 
gle cell suspensions of mLN from infected or immunized mice were rested at 
37°C in a C0 2 incubator for 30 min after isolation. The cells were then stim- 
ulated for 10 min in RPMI-based complete medium in the absence or pres- 
ence of 2.5-10 mM H 2 0 2 for 10 min or with 10 ng/ml OVA-HA for 5 min. 
After stimulation, the cells were immediately fixed by adding a 1 : 1 volume of 
warm Phosflow fix buffer I (BD), followed by incubation at 37°C for 10 min. 
Fixed cells were permeabilized by 30 min methanol treatment on ice, and 
then washed twice with 2% FBS-DPBS and stained with phospho-specific 
antibodies PE— anti-phospho-Erkl/2 (Thr202/Tyr204), PE— a nti— phospho - 
Syk (Tyr525/526), or Alexa Fluor 647-anti-BLNK (pY84) at room tempera- 
ture for 1 h.The cells were then analyzed by flow cytometry. 

Immunohistology. Oct frozen mLNs were processed as described pre- 
viously (Erazo et al., 2007). Sections were stained with antibodies against 
IgE (R1E4), IgGl (RMG1-1), CD35 (8C12), and nuclei were counter- 
stained with Hoechst 33342. Images were acquired with an 1X81 micro- 
scope (Olympus) and Image Pro MDA software (v6.2; Media Cybernetics) 
at 400x magnification. 

Analysis of Sgl remnants in Sjjl-Se junctions. Genomic DNA was ex- 
tracted from sorted IgE + cells and the S(x-Se junctions were amplified by 
PGR using high fidelity ExTaq Polymerase (Takara Bio Inc.) with the prim- 
ers SjaF and SeR (Xiong et al., 2012b). Amplified S|x-Se fragments were 
cloned intoTA vector (Invitrogen), and colonies containing bona fide Su,-Se 
regions were identified by QPCR using primers specific for the (x-Enhancer 
region (MuEnF, 5'-AGCTTGAGTAGTTCTAGTTTCCCC-3' and MuEnR, 
5'-TGGGGACCAATAATCAGAGGG-3') and the Se switch region (SeF, 
5'-GGGCTGCACTAAATGGGACT-3' and SeR, 5'-GCCCGATTGGC- 
TCTACCTAC-3').The presence of switch S7I region footprints was deter- 
mined by QPCR analysis of S7I repeats using S7IF and S7IR primers(Xiong 
et al., 2012b), or by DNA sequencing. 
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VDJ sequence analysis. Total RNA extraction and cDNA synthesis were 
performed using standard procedures with sorted cells isolated from OVA- 
PEP 1— immunized TBmc mice. PGR analysis of cDNA was performed using 
ExTaq Polymerase (Takara Bio Inc.) with the forward primer 5'UTR and 
reverse primers specific for IgE and IgGl constant regions, as previously de- 
scribed (Erazo et al., 2007). 

Adoptive transfer experiments. Switched B cells were isolated from 
pooled spleen and mLN of TBmc mice three weeks after immunization with 
OVA-HA in alum. Cells were preenriched by incubating with negative selec- 
tion cocktail including CD3e-biotin, CDllc-biotin, CD1 lb-biotin, Gr-1- 
biotin, andTER-1 19-biotin, followed by addition of strep tavidin MicroBeads 
and magnetic depletion using the Miltenyi Biotec system. Unbound cells 
were then stained with Alexa Fluor 647 anti-IgE (R1E4), PerCP-eFluro710 
anti-IgGl, PE anti-CD138, APC-eFluro780 anti-B220, PE.Cy7 anti-IgM, 
and eFluro450 anti-IgD. Total switched GC/memory cells were sorted as 
IgM~IgD~CD138~B220 + cells. IgE-switched GC/memory cells were sorted 
as IgE-IgM-IgD _ CD138"B220 + cells. At 3 wk after immunization, IgE 
comprised ^1% of the B220 + switched GC/memory pool. Naive OVA- 
specific CD4 + T from untreated TBmc mice were first enriched by negative 
selection using a CD4 + T cell isolation kit (Miltenyi Biotec) and then stained 
with antibodies against B220 and CD4 for further purification using an 
Aria II cell sorter. Sorted B cells (~2 x 10 6 ) together with 10 5 naive OVA- 
specific CD4T cells were injected i.v.into "y-irradiated BALB/c mice (150 rad 
low dose exposure). The recipient mice, as well control mice that were either 
untransplanted or received T cells only, were subsequently immunized with 
100 jag OVA-HA in alum. IgE-depleted CD138+ plasmablasts/p re-PCs and 
IgE + CD138 + cells were also sorted from immunized TBmc mice and in- 
jected i.v. into -y-irradiated naive BALB/c mice (150 rad exposure). PC re- 
cipient mice were not immunized. Spleen, mLN, BM, and serum were 
collected on day 14 after transfer. To investigate the existence of long-lived 
IgE PCs, serum was collected on day 60 after transfer. Levels of HA-specific 
IgE and IgGl transcripts were analyzed by QPCR using cDNA obtained 
from tissue samples. Primers used for IgE transcript analysis were 17/9DJF1, 
5'-GTACGACGAGAACGGGTTTG-3' and IgE2, and the primers for 
IgGl transcript analysis were 17/9DJF1 and 5 ' - G GAT C C AG AGTTC C A- 
GGTCACT-3' (Erazo et al., 2007). Serum levels of HA-specific antibodies 
were determined by ELISA. 

RNA purification, cDNA synthesis and quantitative PCR analysis 
(QPCR). Tissue RNA was extracted using RNeasy mini kits (QIAGEN). 
cDNA was synthesized with Superscript II reverse transcription (Invitrogen) 
using standard protocols. QPCR analysis was performed as previously de- 
scribed to assess the expression of the Ce switched transcript and Ce and 
C7I mature transcripts for the 17/9 VDJ KI genes ofTBmc mice (Erazo et al., 
2007). The following primers were used for quantification of secreted and 
membrane Ig: secreted IgE, 5'-GTCGCCTAGAGGTCGCCAAG-3' and 
5'-CATCCACCTTCCCCACCACAGC-3'; membrane IgE, 5'-CGCC- 
TAGAGGTCGCCAAGAC-3' and 5 '-AGCTCACACTGAGCAGGAAC-3 ' ; 
secreted IgGl, 5'-TGCACAACCACCATACTGAGA-3' and 5'-GGG- 
TGGAGGTAGGTGTCAGA-3 ' ; membrane IgGl, 5'-TGCACAACCAC- 
C ATAC T GAGA- 3 ' and 5'-CCTCAGCACAGGTCTCGTC-3'. The VDJ 
primers 5'-TGAAACTCTCCTGTGCAGCCTC-3' and 5'-TGCAGAG- 
ACAGTGACCAGAGT-3' were used to quantify total mature IgH transcripts 
in purified IgE + and IgGl + cells from TBmc mice. Expression values in each 
sample were normalized to |3-actin. 

Antigen-specific IgE and IgGl ELISA. HA-specific IgE and IgGl anti- 
bodies were quantified by ELISA as previously described (Curotto de Lafaille 
et al., 2001). To measure HA-specific IgE, serum samples were first de- 
pleted of IgG antibodies by incubation with GammaBmd Plus Sepharose 
(GE Healthcare). 

Mathematical model of GC B cell class switch and implications. An 

agent-based model was developed in which every B cell displayed individ- 
ual dynamics of interaction and development in 3D space and time, with a 



particular emphasis on B cell class switching. The parameters of B andT cell 
migration were set based on data from intravital multi-photon imaging ex- 
periments (Figge et al., 2008). B cell selection in the LZ relies on affinity- 
dependent collection of antigen from follicular DCs and subsequent processing 
and presentation to TFH; hence, affinity-dependent competition forTC help 
is an integral component of the selection process (Meyer-Hermann et al., 
2006;Victora et al., 2010). Zoning of the GC is derived from distinct patterns 
of chemokine distribution {Allen et al., 2004; Meyer-Hermann et al., 2012). 
B cell division and fate decision relies on asymmetric distribution of pre- 
viously collected antigen between the daughter cells (Meyer-Hermann et al., 
2012;Thaunat et al., 2012). All parameters used in these simulations have 
been described in detail elsewhere (Meyer- Hermann et al., 2012). 

The model developed in the current report added a more complex motil- 
ity concept to apoptotic B cells. The phagocytosis rate of apoptotic cells was 
determined by the overall frequency of apoptotic B cells in the experiments. 
The migration model of apoptotic LZ B cells (of any Ig-class) was configured 
to match the experimentally observed DZ/LZ ratio of apoptotic IgGl + cells. 

The switch model used in the current report includes the addition of 
a property of the C++ B cell class. After every division event, class switch 
occurs with a probability given by a transformation matrix (the precise pa- 
rameters are discussed in the main text). Because IgE + GC cells were found 
to express threefold less BCR molecules than IgGl + cells, the probability of 
binding antigen was correspondingly reduced for the IgE + population. Non- 
linear changes of the binding probability were also tested. The lack of BCR 
signaling observed in vivo was implemented in silico by defining the failure 
of IgE + DZ B cells to become sensitized to CXCL13 upon differentiation to 
the LZ phenotype. This assumption was capable of replicating the DZ/LZ 
ratio of IgE + GC cells observed in the experimental data. Each set of simula- 
tion conditions was tested 30 times with different random number seeds and 
then evaluated by calculation of mean and standard deviation. 

Statistics. Statistical significance was determined using Student's t tests 
(unpaired, two-tailed) in Prism 5 (GraphPad Software). 

Online supplemental material. Fig. SI describes the CeGFP KI construct 
and the validation of GFP expression during in vitro induction of CSR to 
IgE. Fig. S2 shows sorting gates and Ig transcript analysis in sorted donor 
IgE + and IgE - cells. A list of antibodies used is provided in Table SI. DZ and 
LZ gene signatures derived from IgGl + GC cells (related to Fig. 6) are pro- 
vided in Table S2. Online supplemental material is available at http://www 
jem.org/cgi/content/full/jem.20131539/DCl. 
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